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We study theoretically the resonant inelastic x-ray scattering in quasi-one-dimensional in- 
sulating copper oxides, where the incident photon energy is tuned to the Culs-4p absorption 
energy. Our attention is focused particularly on the strong momentum-transfer dependence of 
the spectral shape observed in recent experiments. We describe the antiferromagnetic ground 
state within the Hartree-Fock theory, and consider charge excitations from the ground state 
within the random phase approximation. By taking account of the electron correlation effects 
perturbatively, we obtain detailed momentum-transfer dependence of the spectra in a semi- 
quantitative agreement with the experiments. 
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Resonant inelastic x-ray scattering (RIXS) measure- 
ments are now becoming a promising and unique experi- 
mental technique to clarify the detailed charge excitation 
spectra of solids in a relatively high-energy range, 1 owing 
to high brilliance of synchrotron radiation. In particular, 
the RIXS measurements in the hard x-ray regime could 
be a new powerful tool to detect momentum dependent 
charge excitations in solids. A number of compounds 
have been investigated by this technique . 2 ~ 9 Among 
them, the RIXS in cuprates utilizing the Culs-4p absorp- 
tion edge attracts much interest recently, since relatively 
large momentum dependence was indeed observed in the 
charge excitation spectra . 3 ~ 7 > 9 

The electronic properties of the quasi-one-dimensional 
(Q1D) insulating cuprates Sr 2 Cu03 (SC0213) and 
SrCuC>2 (SCO 112) are well characterized by the one- 
dimensional chain consisting of the Cu-0 plaquettes. The 
Cu-0 chains are expected only weakly coupled with each 
other. 10 In SC0213 the Cu-0 plaquettes are aligned with 
sharing the corner oxygen atoms with the nearest plaque- 
ttes, as shown in Fig. 1. It was reported that SC0213 
shows antiferromagnetic order below about 5 K with 
reduced staggered magnetization along the chain. The 
strong fluctuations are responsible for the extremely re- 
duced Neel temperature and magnetization, 11 but the 
ground state of SC0213 is an antiferromagnetic (AF) in- 
sulator. Concerning SC0112, the Cu-0 plaquettes form 
a zigzag chain by sharing the edges with the nearest pla- 
quettes. The zigzag chain can be regarded approximately 
as being constructed by combining two independent Cu- 
O chains sharing the edges with each other. 10 Therefore 
we study the RIXS properties of these two cuprates by 
considering the corner-sharing Cu-0 chain as in Fig. 1 in 
the present article. 

Recently, RIXS measurements were performed for the 
Q1D cuprates SC0213 6 and SC0112. 6 ' 9 In the both Q1D 
cuprates, characteristic spectral weight was obtained for 
energy transfer u> » 2,5.6 (eV), but the origin of the 
spectral peaks are still unclear. The aim of the present 
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Fig. 1. Schematic figure of the Cu-O chain presenting the defini- 
tion of the transfer integrals. 



article is to provide a quantitative understanding on 
the RIXS, particularly on the position of the peaks and 
the highly momentum-dependent behaviors of the spec- 
tra in these one-dimensional cuprates. We calculate the 
intensity in the RIXS spectra as a function of energy- 
momentum transfer, using the perturbation theory de- 
veloped by Platzman and Isaacs. 14 We describe the an- 
tiferromagnetic (AF) ground state by the Hartree-Fock 
(HF) theory, and take account of the excitations from 
the ground state within the random phase approxima- 
tion (RPA). As a result, we obtain a semiquantitative 
agreement with experiments, regarding the peak posi- 
tions and the momentum dependent behaviors. 

The essential microscopic process of the RIXS in the 
cuprates is considered to be as follows. 12 The incident 
photon excites resonantly the Culs core electrons to 
the upper empty Cu4p bands. Since the core Culs or- 
bitals are well localized to be coupled strongly to the 
Cu3d electrons by the interorbital Coulomb interaction, 
the created Culs core hole scatters the Cu3d electrons, 
which occupy the bands relatively close to the chemi- 
cal potential (As a result of the scattering, the most 
probably only one electron-hole pair is created on the 
Cu3<i band). In order to describe this process specif- 
ically, we consider the total Hamiltonian of the form, 
H = Hdp-\-H ls _3d+H ls +H± p + H x . The dp-Hamiltonian 
Hd p describing the electronic properties of the Cu3d-02p 
bands is given in the form, Hd p = Hq + H' . Hq and H' 
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are the nonintcracting and the interaction parts, respec- 
tively. 
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where d ka = {d ka ,Pika,P2ka,P2'ka), and d ka , p lka , p 2kiy 
and P2'fccr are the electron annihilation operators in the 
atomic d x 2_ y 2 or p x y orbits at the Cu, 0(1), 0(2) and 
0(2') atoms, respectively. The dispersions are CiCO — 
2itisin-| and C'(^) = 2ii'sin-|. The hopping parameters 
are h = -1.45 eV, t 2 = -1.8 eV and t' = -0.7 eV 
(See Fig. 1), as determined for SC0213 by the LDA first 
principle calculation. 13 The interacting part is 

tf' = ^£4- 9 t4' +g Ai4 T . (2) 

kk' q 

The Coulomb energy is taken to be U,j = 11 eV in the 
present study. Regarding the one-particle level, we take 
s P2 — s pi = 0.5 cV, e^ F — e pi = —0.5 eV, where ejp 
is the screened Cu3<i one-electron energy within the HF 
theory. For this parameters, we obtain the AF ground 

state with the staggered magnetization m s tagg. = 0.43/ie- 
The scattering of the Cu3d electrons by the Culs core 
hole is described by 

H u-3d = £ d ] k , +qa s\_ qa ,s k(T ,d kl(T , (3) 

kk' qua' 

where s ka (sjL-) 1S ^ nc annihilation (creation) operator 
for the Culs electrons with momentum k and spin a, 
and V is the so-called core-hole potential. H^ v and H\ s 
describe the kinetics of the electrons on the Cu4p and 
Culs bands. Since the Culs electrons are well localized, 
we take completely flat dispersion for them. For the Cu4p 
electrons, we use simple two-dimensional cosine-shaped 
band for simplicity. This simplification does not affect 
the spectral shape so drastically, because the factor con- 
taining the Cu4p dispersion function is integrated up in 
momentum, as seen later in eq. (5). H x describes the 
transitions between the Culs and Cu4p states, involving 
the photon absorption and inverse emission processes. 
H x is of the form 

H x =^2(w(q;k)p' k \ qa s ka +h.c), (4) 

kqcr 

where p'^ is the creation operator of the Cu4p electron. In 
the present study, we ignore the momentum dependence 
of the matrix elements w(q; k), i.e., w(q; k) = w. 

The main spectral weight of the RIXS is evaluated 
by the diagram shown in Fig. 2. In Fig. 2(a), the effec- 
tive three point vertex function A(w,q) represented by 
the shaded triangle is calculated within the RPA, and 
this vertex part is inserted to the diagram Fig. 2(b) rep- 
resenting the total scattering probability. In Fig. 2(b), 
note that the off-diagonal components of the Keldysh 
Green's functions are assigned to the solid lines which 
connect the upper normally-time-ordered branch and 
the lower reversely-time-ordered branch, while the usual 
causal (normally-time-ordered) Green's functions are as- 
signed to the solid lines in the upper branch and the 



(a) 




Fig. 2. (a) Vertex renormalization in the RPA. The solid lines, 
the shaded triangle and the empty square are the Green's func- 
tion for Cu3d electrons, the renormalized vertex and the anti- 
symmetrized bare Coulomb interaction, respectively, (b) The di- 
agrammatic representation of the RIXS intensity. The wavy and 
the dashed lines denote photons and the ls-core-holc potential, 
respectively. The Green's functions for Cu3d, Culs and Cu4p 
electrons are assigned to the solid lines marked with '3d', 'Is' 
and '4p', respectively. 



reversely-time-ordered Green's functions are assigned to 
the solid lines in the lower branch. 15 Regarding the core 
hole potential V, we take only the first order term in V 
(i.e., Born scattering). The effect of higher orders in V is 
briefly mentioned later. The analytic expression for the 
diagram Fig. 2(b) is obtained as 

Wiwm-^fQf) = (2x) 3 AT|«,| 4 ^ S(E j {k)+u-E j ,(k + q)) 

kjj' 

xn 3 (fc)(l - n jt (k + 9))hr U^aA^K*' W, 9)C / I„/,y (* + l) 



4^ (u>» + Els + iri s - e4p(fel))(wf + Els + iTl B — E4p(fel)) I 

where qi and qf (oji and ojf) are the momenta (energies) 
of the initially absorbed and the finally emitted photons, 
respectively, q = qi — qf and lo = Ui — ujf are the momen- 
tum transfer and the energy loss, respectively. Ej(k) and 
rij(k) are the energy dispersion and the electron occupa- 
tion number of the band j, respectively, obtained by di- 
agonalizing the HF Hamiltonian for Hd p . Uj^da{k) is the 
(J, (icr)-element of the unitary matrix diagonalizing the 
HF Hamiltonian for Hd p - £i s and £4 P (/c) are the kinetic 
energies of the Culs and Cu4p electrons, respectively. In 
the present study the incident photon energy u>i is tuned 
to the Culs-4p absorption energy Wj w £4 P (0) — e\ s , for 
which the intensity of the spectrum is enhanced by the 
resonance. The decay rate Ti s of the core hole is 0.8 eV 
in the present study. 

The numerical results of the RIXS intensity by the 
formula eq. (5) are shown as a function of energy loss 
lo for momentum transfers q — 0, § and 7r in Fig. 3. 
We see that there are three characteristic peaks in the 
calculated RIXS spectra around lo w 2, 6, 10 eV. In the 
experiments, two peaks are observed around the momen- 
tum transfer lo w 2 and 5.6 eV, which are in agreement 
with the present numerical results. The high-energy lOeV 



J. Phys. Soc. Jpn. 



Letter 



Takuji NOMURA and Jun-ichi IGARASHI 



3 



• iH 

c 

1 



C 

1) 



C/3 



q=0 





q=7U/2 



q=7t 




4 6 8 10 

Energy Loss (eV) 



12 



Fig. 3. The calculated spectra as a function of energy loss for 
three momentum transfers q = 0, 0.5-7T and 7r. The thick solid and 
the thin broken lines are the results with and without the RPA 
corrections, respectively. The inset shows a typical experimental 
spectrum read from ref. 6, and the arrows point the two peaks. 



peak gives only small contribution to the total spectral 
weight, and the intensity around uj « 2 and 6 eV cov- 
ers the main part of the spectral weight. In Fig. 4, the 
contourplots of the calculated RIXS spectra are shown 
as a function of the momentum transfer q and the en- 
ergy loss to, and the experimental data read from ref. 6 
are also displayed for comparison. We can see that the 
2eV peak shows a relatively large sinusoidal dispersion, 
while the other two do not. This sinusoidal dispersion 
of the 2eV peak quantitatively agrees with the experi- 
ments. 6,9 Concerning the 6eV peak, the calculated peak 
positions deviate somewhat from the experimental ones, 
maybe due to the crude tight-binding fitting. The weak 
momentum dependence of the 6eV peak position is well 
reproduced by the calculation. 

In order to inspect the origin of the peaks in the RIXS 
spectra, we should turn back to the formula (5). Com- 
paring the results modified by the correlations (the solid 
lines in Fig. 3) with the uncorrected results (the dashed 
lines in Fig. 3) , we can see that the overall spectral shape 
is roughly captured even without any correlation effects, 




71/2 71 37C/2 

Momentum Transfer q 



2tu 




tc/2 

Momentum Transfer q 



Fig. 4. Graylevel plots for the spectra as a function of momen- 
tum transfer and energy loss. The light and dark regions corre- 
spond to high and low RIXS intensity, respectively. The hatched 
squares and circles are experimental peak positions in SC0213 
and SCOH2, respectively, from ref. 6. In (b), the diamonds con- 
nected by a thin line denote the theoretical peak positions. 



although the detailed spectral structure are of course af- 
fected by the correlation effects. Therefore we might ex- 
pect that the overall spectral shape is determined by the 
other factors rather than the correlation effects on the 
scattering process, and thus expect naturally that the 
factor rijik)^ - rij, (k + q))\Uj^(k)U f da , d ,(k + q)\ 2 is es- 
sential for determining the overall spectral shape. It is 
very interesting to note that this factor is the product 
of the partial electron occupation number of the band 
j at the momentum k [ nj(k)\Uj,da(k)\ 2 ] and the par- 
tial hole occupation number of the band j' with the mo- 
mentum k + q [ (1 - rij>{k + q))\U 1 d<T , Jt (k + q)\ 2 ]. Thus 
the Culs-4p RIXS spectra in cuprates are closely related 
to the partial occupation number of the Cu3d electrons 
in each band. From this point of view, the RIXS spec- 
tra reflect the charge excitation processes in which Cu3e? 
electrons are selectively excited from the lower occupied 
bands to the upper unoccupied bands. In Fig. 5, the total 
and the Cu3ci partial density of states (DOS) are shown. 
We find only one DOS weight around to ~ I eV above 
the chemical potential, and four DOS weights around 
to « —2, —4, —6, —9 eV below the chemical potential. 
We can consider naturally that the three RIXS weights 
around 2, 6 and 10 eV in Fig. 3 are attributable to the 
excitation processes in which the electrons occupying the 
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Fig. 5. The density of states for the antiferromagnetic ground 
state. The thin dotted and thick solid lines represent the total 
and the Cu3d partial density of states, respectively. The chemical 
potential is set to ui = 0. The inset shows the band dispersions. 



bands u> « —2, —6, and —9 eV are excited to the upper 
unoccupied band. Of interest is that the sharp peak in 
the total DOS around w ~ — 4 eV, which corresponds 
to the non-bonding oxygen band, is not reflected in the 
RIXS spectra at all. This is because the oxygen non- 
bonding band does not contain any components of the 
Cu3d orbitals and the other orbitals which are enough 
localized at the Cu-sites to be coupled strongly to the 
Culs orbitals. In this sense, the Culs-4p RIXS measures 
the Cu3d-orbital-selective charge excitations. In general, 
only the bands which are hybridized with the localized 
orbitals strongly coupled to the core hole are observable 
in the RIXS spectra. 

We should make a distinction between the 2, 6 eV 
peaks and the 10 eV peak in the RIXS spectra (Fig. 3). 
Since the DOS peaks around ut w 1 and —9 eV cor- 
respond respectively to the upper and lower Hubbard 
bands, the RIXS spectral weight around 10 eV in Fig. 3 is 
related to the charge excitation process across the Mott- 
Hubbard energy gap, which is usually of the order of Ud- 
On the other hand, the RIXS spectral weights around 2 
and 6 eV are related to the charge transfer (CT) pro- 
cess in which the Cu3(i electrons mixed in the wide 02p 
bands are excited to the upper empty Cu3d Hubbard 
band. This energy scale of the CT excitation is basically 
connected to the CT energy e p — e<j. We should note 
that the main spectral weight originates from the CT 
process, and the excitations across the Mott-Hubbard 
gap provides only minor contribution to the RIXS spec- 
tra. Therefore, the validity of the theoretical analyses 
based on the single-band Hubbard model, 12 which con- 
sider only Cu sites, is unclear. 

We give some remarks on the present calculations. At 
first, we discuss the validity of the application of the HF 
theory to the present system. Although the experiments 
are usually performed above the Neel temperature, the 
HF analysis is valid in principle only at the absolute zero 
temperature. In actual low-dimensional systems includ- 
ing the cuprates, the fluctuations are very strong to re- 
duce the Neel temperature and the magnitudes of order 
parameters. However we would like to stress that it is 



not so important for the RIXS spectra whether the AF 
long range order is attained or only short range AF cor- 
relation exists. The reasons are as follows, (i) The AF 
transition does not affect the RIXS spectra so drasti- 
cally, since the RIXS spectra basically reflect only the 
charge sector in the scattering process, (ii) The RIXS 
occurs, basically being well localized in the space, due to 
the localized nature of the Culs core hole. In addition 
to these points, by noting that the AF ground states in 
Mott insulators may be continued adiabatically to the 
HF AF ground state at the absolute zero temperature 
in the limit of strong Coulomb repulsion, we consider 
that the AF HF theory provides a good starting point 
to the analysis. 16 But we should note here the following 
point. As we have mentioned, the RIXS spectral weight 
is roughly determined by the product of the Cu3rf partial 
occupation numbers in the occupied and the unoccupied 
bands. Since the HF theory usually underestimates the 
magnitude of the localized moment, the product of Cu3d 
electron and hole occupation numbers rii a (l — rtfo) in the 
local picture is overestimated. Thus the intensity around 
10 eV in the RIXS spectra is overestimated within the 
HF theory. It is natural to consider that the CT exci- 
tations, rather than the charge excitation process across 
the Mott-Hubbard gap, provide the main contribution 
to the spectral weight, and the 10 eV peak might not be 
observable in reality. 

We have considered only the first order terms in the 
core hole potential V. Actually we investigated the ef- 
fects of higher orders in V within the T-matrix approxi- 
mation as in usual single impurity problems (Results are 
not shown). However we found that the spectral shape is 
not so drastically changed by the higher orders of V, al- 
though the absolute magnitude of the spectra is changed. 

In the present work, we have analyzed theoretically 
the RIXS spectra in the Q1D cuprates. The peak po- 
sitions and the momentum-transfer dependence of the 
RIXS spectra are explained in a semiquantitative agree- 
ment with the experimental results. We have shown that 
the RIXS measurement in cuprates is a unique technique 
to clarify the Cu3(i-orbital-selective charge excitations. 

The authors are much grateful to Prof. M. Hasan for 
the communications. The numerical computation was 
partly performed at the Yukawa Institute Computer Fa- 
cility. 
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